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Abstract 

In the coevolving voter model, each voter has one of two diametrically opposite opinions, and 
a voter encountering a neighbor with the opposite opinion may either adopt it or rewire the 
connection to another randomly chosen voter sharing the same opinion. As we smoothly change 
the relative frequency of rewiring compared to that of adoption, there occurs a phase transition 
between an active phase and a frozen phase. By performing extensive Monte Carlo calculations, we 
show that the phase transition is characterized by critical exponents (3 = 0.54(1) and v = 1.5(1), 
which differ from the existing mean-field-typed prediction. We furthermore extend the model by 
introducing a contrarian type that tries to have neighbors with the opposite opinion, and show that 
the critical behavior still belongs to the same universality class irrespective of such contrarians' 
fraction. 
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I. INTRODUCTION 



Analyzing social dynamics from a physical viewpoint has become an important topic in 
statistical physics. In particular, emergence of collective opinion is an intriguing phenomenon 
which can be readily compared to the collective ordering manifested in many statistical- 
physical systems, and the voter model has served as an insightful starting point to study 
the opinion dynamics [l| . When the voter model began to be investigated in the context of 
complex networks, the underlying connection structure among voters was usually assumed 
to be static. In this case, the dynamics can be basically described as follows: The opinion 
Sj assigned to a voter j is either +1 or —1 like an Ising spin. For a given voter j, we 
randomly choose one of the neighbors, say k. If Sj ^ Sk, the link between j and k is 
called active, and j adopts /c's opinion under certain stochastic rules. As implied in this 
adoption mechanism, it is true that our opinion is shaped by whom we are surrounded by, 
but it can be equally true that our opinion shapes our social networks. So the situation 
becomes more interesting when we allow voters to change their neighbors. It then looks 
natural to introduce a tendency to connect to like-minded people, because as the proverb 
says, birds of a feather flock together. This simultaneous change in the network structure is 
termed coevolution, and the voter model has been studied in this coevolutionary framework 
for the natural reason mentioned above [2t5|. Besides the opinion dynamics, we can also 
mention that the coevolutionary dynamics plays a crucial role in the context of epidemiology 
since adaptation of a network structure in response to disease spreading is directly related 
to quarantine policies (see Refs. for theoretical investigations). In addition, more 



physically motivated studies include a coevolving network of phase oscillators 
Ising model with Glauber-like dynamics jol , and a variation of the XY model 



the kinetic 



io|. 



II. COEVOLVING NETWORK OF VOTERS 

In this work, we begin with the coevolving voter model proposed by Vazquez et al. ^ 
where an active link becomes inert either by rewiring with probability p or by adoption with 



probability 1 — p (see also Ref. [ll| for its extended versions). According to a mean- field 
(MF) argument and numerical simulations, this model has a well-defined phase transition at 
a critical value Pc, separating an active phase from a frozen phase in which one observes no 
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FIG. 1. (Color online) Possible states of a link between a pair of voters. We represent a conformist 
as a filled triangle, pointing up (a) or down (t) depending on the opinion. Likewise, a contrarian 
is represented as an empty triangle (A or V). A full red arrow and a dashed blue arrow mean an 
active link and an inert one, respectively. Whether a link is active depends on whose viewpoint is 
taken, so an arrow from a voter means the state of the link from her viewpoint. For example, the 
link between V and A is active from A's viewpoint, while it is inert from V's viewpoint (bottom 
left). 



active links between the two different opinions. Specifically, as p approaches the critical point 
Pc from below, the active-link density is expected to scale as p ~ {Pc~ pY ■, equivalently, 
p ~ N'^/" at p = pc for a large but finite number of voters, where (3 and u are critical 
indices. The MP argument in Ref. ^ predicts /3 = 1 and pc = (/i — 2)/(/i — 1), where n 
means the average degree. Although this predicted significantly deviates from numerical 
estimates, there are better methods to estimate Pc 12|. However, it has not been questioned 
whether the critical behavior is consistent with the MF prediction, and we will show that 
this is not the case by finding /3 = 0.54(1) and u = 1.5(1). Then, we extend the model by 
introducing a peculiar voter type tha t tries to surround herself with the opposite opinion, 

. The effect of such behavior has already been studied 



which might be called heterophily 13 

n n n 

in opinion dynamics [1^ 115| and Ref. [16[ 



has argued that such behavior can be generated by 



stochasticity. Though this heterophilic dynamics is sometimes dubbed the "inverse voter" 
model l^j, let us denote the usual voter and the inverse voter as a conformist (v) and a 
contrarian (v), respectively, following the terminology in Refs. fl^. Iisl. Our point is that 
the critical behavior is universal regardless of such contrarians' fraction whereas Pc is a 
nonuniversal quantity. 

For simplicity, we assume that voters do not change their types, (i.e., once a contrarian, 
forever a contrarian) so the fraction r of conformists is constant in time. When r = 1, our 

n 

model is reduced to the original one in Ref. [2[, while for r = 0, it is identical to the inverse 



voter model in Ref. [17|. It is important to note that whether a hnk is active between two 
voters depends on their types: Fig. [1] shows all the possible link configurations between two 
voters, where the signs mean their opinions. Although each link is expressed as two directed 
arrows, there exists only one single link between a pair of voters so the double arrows imply 
that the state of a link is not uniquely defined, in general. For example, if a conformist with 
opinion +1 is linked to a contrarian with —1, the conformist regards the link as active while 
the contrarian regards it as inert. If a link connecting two voters is regarded by both as 
active (inert), the link contributes 1 (0) to the total number of active links in the system. 
If only one of them regards it as active, on the other hand, it is counted as one-half. 

Let us now explain details of our calculation. In terms of a network, a voter is mapped 
to a node, and a link connecting nodes j and k can be denoted as {j,k). As an initial 
condition, we construct a regular random network of size where every node has exactly 
/i neighbors (we choose /i = 4 throughout this work). We then randomly select rN nodes 
as conformists while the others are set to be contrarians. Every node has equally random 
probability to have either opinion at the starting point. Then, we randomly pick up a node 
j and choose one of its neighbors, say k, also at random. If the link (j, k) is inert from j's 
viewpoint, nothing happens. If it is not, there are two possible options: With probability p, 
[j, k) is rewired to a randomly chosen node / in such a way that {j, I) becomes inert from j's 
viewpoint. Or with probability 1 — p, Sj is fiipped to make (j, k) inert from j's viewpoint. 
At each time step t, this updating procedure is repeated for A^ times. 



III. RESULTS 

We take the active-link density p as an order parameter whose dependence on N at p = pc 
provides estimates of the critical indices. In order to find its stationary value for each A^, 
p should be first measured as a function of t over A^, independent samples. Note that we 
measure surviving average by including only samples with at least one active link at each 
given time t. This surviving-averaged quantity shows transient behavior at small t, whereas 
it fiuctuates violently at large t as statistics become poor with few surviving samples. So the 
stationary value p{N) is estimated by taking an average over a certain time interval [ti,t2]- 
In our case, the lower bound ti is chosen as 5 x 10^ or 10'^ depending on A^, while ^2 is fixed 
at 5 X 10'^. The error bar becomes generally smaller as A^g increases, and we find it quite 
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FIG. 2. (Color online) Coevolving voter model with r = 1 where all voters are conformists, (a) 
The stationary active-link density p is plotted as a function of at different rewiring probabilities. 
The solid line is a guide to eyes, (b) The function / in Eq. ([1]) is obtained by scaling collapse with 
= 1.5(1), /3 = 0.54(1), and Pc = 0.383(1). 

enough to work with A^j, = 10^ for r = 1/2 and A^^, = 2 x 10^ for r = or r = 1. In addition, 
as for the system size A^, we need to consider two limitations: That is, if A^ is too small, it 
is difficult to generate many surviving samples at large t since p = is an absorbing state 
which can always be reached due to fluctuations in a finite-sized system. If A^ is too large, 
on the other hand, it is difficult to approach a stationary value, especially near the critical 
point, due to the diverging relaxation time scale. For these reasons, we will concentrate on 
precise estimates of p{N) from A^ = 10"^ to A^ = 2 x 10'^. 

We start with r = 1 and plot numerical results of p{N) at various values of p in Fig. [2](a). 
It is observed that p decays algebraically at p ~ 0.383, indicating the existence of a con- 
tinuous phase transition at pc = 0.383(1), in agreement with pc ~ 0.38 found by numerical 
simulations in Ref . ^ . From the standard finite-size scaling form 

p = N-^'^f{eN^/^) (1) 

with e = {p —PcilPc-, the slope of the line in Fig. [2|^a) determines /3/z/ 0.36. We use 
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FIG. 3. (Color online) (a) Scaling collapse for r = where all voters are contrarians and (b) for 
r = 1/2. The critical exponents /3/z/ = 0.37(2) and z/ = 1.5(1) are used in common, whereas 
Pc = 0.383(1) in panel (a) in contrast to pc = 0.790(1) in panel (b). 

these values of Pc and /3/z/ to make the data collapse in Fig. [21(b). This analysis finds that 
Pc = 0.383(1), u = 1.5(1), and (3 = 0.54(1). It is noteworthy that the actual value of (3 is 
far from the MF prediction that /3 = 1. The reason is that the full mixing assumption of 
the MF approach cannot be justified when the system splits into clusters, which is exactly 
what happens at the critical point [see Fig. |5]^a) below]. In this respect, it should not be 
surprising that the actual transition is not of the MF type. 

We then repeat the same analysis for r = [Fig. [3t^a)]. The resulting exponents /3/z/ = 
0.37(2) and u = 1.5(1) as well as Pc = 0.383(1) are in striking agreement with those in the 
above case with r = 1. From the observation that both with r = and r = 1 belong to the 
same universality class characterized by the same critical exponents /3 and u, one might guess 
that this holds true for any value of r G [0, 1]. In order to check this, one could perform the 
same finite-size scaling analysis for r = 1/2 [Fig. [3]^b)]. Unfortunately, finite-size corrections 
to scaling are so substantial in our data that we should be content with a consistency check 
just by tuning p^ with having fixed /3/z/ = 0.37 and u = 1.5. In spite of such limitations, we 
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FIG. 4. Phase diagram of the model on the (p, r) plane. See text for details. 

can say that the result is not inconsistent with the guess. It is somewhat interesting to note 
that the coevolving model in Ref. j4| has 1/z/ 0.7(1), which is also close to our finding 
Ijv ^ 0.67. 

Let us now turn our attention to determining the critical threshold to find the shape of 
the phase boundary. The estimate of Pc in the simple MF approximation by Vazquez et 
al. |2| can be briefly explained in the following way: Suppose that there is one active link in 
the system with average degree /i. For example, there are two clusters of opposite opinions 
and there is one single link between them. If the link is rewired with probability this 
active link will disappear. If one voter changes the opinion with probability 1 — p, then the 
original active link disappears but all the other /i — 1 links become active. Therefore, the 
expected change in the number of active links is [{fi — 1) — 1](1 —p) —p, which is zero below 
the threshold, leading to p^ = (/i — 2)/(/i — 1). Note that the above simple MF approximation 
does not depend on the value of r, and thus giving = 2/3 for /i = 4 irrespective of r, 
which is sharply against our observation that Pc{r = 1/2) 7^ pc{r = 1). 

In order to numerically obtain the actual phase diagram on the {p, r) plane, we first 
measure p* = p{p = Pc = 0.383) at r = 1 for size = 10'^ and then use the same p* to 
locate Pc for other values of r. In other words, Pc{r) is roughly estimated from the condition 
that p{pc', r,N = 10^) = p* . We have confirmed that this simple method yields p^ reasonably 
close to the the above-mentioned estimates for r = and r = 1/2. Furthermore, it is enough 
to use a smaller ensemble (A^^ = 200) for the present purpose. The resulting phase diagram 
shows a lobelike shape with a maximum value around r = 1/2 (Fig. Hj), and the MF theory 
fails to predict this feature once again. 

In order to understand such a shape, let us look into the MF estimate a little further: 



In updating a node with the probabihty 1 — it is possible to come back to the original 
configuration by flipping the same node once again. This possibility is ignored in the simple 
3], which means that the spatial correlation among the newly created 



MF scheme in Ref. 



active links is ignored. The g-fan motif method [1^ is a more advanced approximation in the 
sense that it evolves the configuration until clusters of the active links are separated by at 
least one spacing. If the configuration is expressed with a finite number of separate clusters 
of active links after some time, then the method ignores the spatial correlation among them 
in describing their further evolution with a stochastic matrix. Comparing r = with r = 1, 
we see that the g-fan motif method works in the same way and the resulting estimates of Pc 
are therefore identical again. When < r < 1, however, there arises an important difference 
that there can exist an active link between a conformist and a contrarian, which never 
disappears just by changing one's opinion. There is always a finite chance that a cluster of 
active links keeps growing from this active seed link by fiipping neighboring voters. Thus 
the configuration cannot be generally a collection of finite separated clusters for < r < 1, 
and the g-fan motif method fails in estimating pc- Nevertheless, it tells us that one needs a 
higher p in order to suppress such growth of the cluster by rewiring the active link instead 
of changing opinions. This qualitatively explains why the phase boundary at < r < 1 is 
pushed to higher p in Fig. |H 



Finally, we visualize actual network structures in Fig. El When everyone is a conformist 
(r = 1), one can clearly see that the network structure radically changes as p passes the 
critical point so that the frozen phase is characterized by completely segregated opinions. 
More specifically, as Pc is approached from below, two clusters have fewer and fewer inter- 
cluster links and they get disconnected at pc- Therefore, in the vicinity of the critical point, 
the full-mixing assumption in the simple MF approximation just breaks down, allowing the 
critical exponents to deviate from the MF ones, as confirmed in the present work. We 
add a remark that network structures in the frozen phase largely depend on the conformist 
fraction r. For example, when r = 1/2, the system usually splits into three clusters, one for 
contrarians and the other two for conformists. 
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(a)r=l,p=0.3 {b)r=l,p=0.5 




FIG. 5. (Color online) Snapshots of network structures in the active phase [(a), (c), (e)] and 
in the frozen phase [(b), (d), (f)] with r = 1 [(a),(b)], r = 1/2 [(c), (d)], and r = [(e), (f)]. The 
number of nodes is = 200 in each case. Circles and triangles denote conformists and contrarians, 
respectively, and a filled (empty) symbol denotes the opinion +1 (—1). Red, blue, and gray colors 
mean that the link is active in both directions, active in only one direction, and inert in both 
directions, respectively. All links become inert (gray) in the frozen phase. 



IV. SUMMARY 



In summary, we have numerically studied the phase transition of the coevolving voter 
model composed of conformists and contrarians. Our main result suggests that the simple 
MF approximation cannot explain the transition nature, and that the universality class is 
well characterized by critical exponents f3 ^ 0.54 and u ~ 1.5, irrespective of the conformist 
fraction. We present the full phase diagram on the (p, r) plane, which shows a symmetric 
lobelike structure around r = 1/2. 

9 



ACKNOWLEDGMENTS 



B. J.K. was supported by the National Research Foundation of Korea (NRF) grant funded 
by the Korea government (MEST) (Grant No. 2011-0015731) and C.-RZ. by the Natural 
Science Foundation of China (NSFC) (Grant No. 11175086). 



[1] See, e.g., C. Castellano, S. Fortunate, and V. Loreto, Rev. Mod. Phys. 81, 591 (2009) for 
further references. 

[2] F. Vazquez, V. M. Eguiluz, and M. San Miguel, Phys. Rev. Lett. 100, 108702 (2008). 

[3] R. R. Durrett, J. P. J. Gleeson, A. L. A. Lloyd, P. J. P. Mucha, F. F. Shi, D. D. Sivakoff, 

J. E. S. J. Socolar, and C. C. Varghese, Proc. Natl. Acad. Sci. U.S.A. 109, 3682 (2012); 

I. Benczik, S. Benczik, B. Schmittmann, and R. Zia, Phys. Rev. E 79, 046104 (2009); F. Fu 

and L. Wang, ibid. 78, 016104 (2008). 
[4] P. Holme and M. E. J. Newman, Phys. Rev. E 74, 056108 (2006). 
[5] G. Ihiguez, J. Kertesz, K. Kaski, and R. Barrio, Phys. Rev. E 83, 016111 (2011). 



[6] T. Gross, C. J. D. D'Lima, and B. Blasius, Phys. Rev. Lett. 96, 208701 (2006) T. Gross and 



B. Blasius, Journal of The Royal Society Interface 5, 259 (2008) N. Crokidakis and S. M. D. 



Queirs, |J. Stat. Mech: Theory Exp. 2012, P06003 (2012) 



[7] L. B. Shaw and I. B. Schwartz, Phys. Rev. E 77, 066101 (2008) 



[8] T. Aoki and T. Aoyagi, Phys. Rev. Lett. 102, 034101 (2009). 



[9] S. Mandra, S. Fortunato, and C. CasteUano, Phys. Rev. E 80, 056105 (2009) K. B. Hajra 



and A. K. Chandra, Eur. Phys. J. B 85, 27 (2012). 



[10] P. Holme, Z.-X. Wu, and P. Minnhagen, Phys. Rev. E 80, 036120 (2009) 



[11] J. Cheng, Y. Hu, Z. Di, and Y. Fan, Comp. Phys. Comm. 181, 1697 (2010); J. L. Herrera, 
M. G. Cosenza, K. Tucci, and J. C. Gonzalez-Avella, FPL 95, 58006 (2011); C. Nardini, 
B. Kozma, and A. Barrat, Phys. Rev. Lett. 100, 158701 (2008); G. Zschaler, G. A. Bohme, 
M. Seifiinger, C. Huepe, and T. Gross, Phys. Rev. E 85, 046107 (2012). 

[12] G. A. Bohme and T. Gross, Phys. Rev. E 83, 035101 (2011). 

[13] D. Kimura and Y. Hayakawa, Phys. Rev. E 78, 016103 (2008). 

[14] S. Galam, Physica A 333, 453 (2004). 

10 



[15] N. Crokidakis and C. Anteneodo, |Phys. Rev. E 86, 061127 (2012)t S. Biswas, A. Chatterjee, 



and P. Sen, [Physica A 3 91, 3257 (20 12) 



[16] M. S. De La Lama, J. M. Lopez, and H. S. Wio, EPL 72, 851 (2007). 

[17] C.-P. Zhu, H. Kong, L. Li, Z.-M. Gu, and S.-J. Xiong, Phys. Lett. A375, 1378 (2011). 

[18] H. Hong and S. Strogatz, Phys. Rev. Lett. 106, 054102 (2011). 



11 



